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We carry out a detailed examination of the ground state property of few-boson system in a 
one-dimensional hard wall potential with a ^— split in the center. In the Tonks- Girardeau limit 
with infinite repulsion between particles, we use the Bose-Fermi mapping to construct the exact 
A/"— particle ground state wavefunction which allows us to study the correlation properties accu- 
rately. For the general case with finite inter-particle interaction, the exact diagonalization method 
is exploited to study the ground-state density distribution, occupation number distribution, and 
momentum distribution for variable interaction strengths and barrier heights. The secondary peaks 
in the momentum distribution reveal the interference between particles on the two sides of the split, 
which is more prominent for large barrier strength and small interaction strength. 

PACS numbers: 03.75.Mn, 03.75.Hh, 67.60. Be 



I. INTRODUCTION 

With the development of laser cooling and optical trap- 
ping technology, quasi one-dimensional (ID) cold atom 
systems have been realized by tightly confining the par- 
ticle motion in two directions to zero point oscillation 
P, in, S 0, Q- Furthermore, optical dipole forces gen- 
erated by several crossed, interfering laser beams make it 
possible to create periodic potentials, and experimental- 
ists are able to trap and control small numbers of parti- 
cles in optical lattice or double well potential. Meanwhile 
with the Feshbach resonance technique, the inter-particle 
scattering length can be easily changed by tuning a mag- 
netic field, which allows the atoms to enter the full regime 
of interaction from weakly to strongly interacting limit. 
In particular, a quasi- ID quantum gas of strongly inter- 
acting bosons has been observed in the so-called Tonks- 
Girardeau (TG) regime [3|, y. These experiment pro- 
gresses inspire great interest in exploring the properties 
of bosonic gases in the full interaction regime. 

As is well known, the Gross-Pitaviskii (GP) theory 
is widely adopted in dealing with the system of Bose- 
Einstein condensations with weak interaction. In spite 
of its great success in explaining the basic experimental 
observations, the GP equation is fundamentally based 
on a mean-field approximation and fails in the strongly 
correlated systems 0, Q • For arbitrary interaction from 
weakly to strongly interacting limit, there exist sev- 
eral theoretical treatments including analytic method, 
e.g. Bethe ansatz [ol, [Tol, [O, EH, numerical 
simulation e.g. the exact diagonalization [3, [lH, the 
multi-configuration time-dependent Hartree (MCTDH) 
method OIl^? discrete variable representation (DVR) 
[l^, and so on. Bethe ansatz provides the exact so- 
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lution to some many- body systems including Bose gas 
Fermi gas [ll|, ll2[, Bose-Fermi mixture [13] as 
well as multi-component cases. In the past several years, 
though it has made great development in cooperation 
with other theories and some numerical methods, this 
analytical method mainly solves the uniform system. On 
the other hand, numerical techniques provide more in- 
formation about the many-body system, for example, 
MCTDH method can not only reveal the ground state 
property, but also tackle its dynamics [19]. For finite 
particle interactions no analytical solution to the two- 
particle problem is known for a finite split barrier in the 
center. DVR has been used to achieve the discretization 
of the spatial coordinates for two particles, however, the 
generalization to more particles are not straightforward 

While most theoretical works assume large particle 
numbers, strongly correlated atomic systems realized in 
current experiments are limited to systems with a small 
number of particles 0, 0]. The exact diagonalization 
method becomes more and more versatile and powerful 
in dealing with the ground state problem of such "few- 
body" system. In this paper, we study the ground state 
properties of a few spinless bosons which are confined in 
a ID hard- wall trap and interact via repulsive contact po- 
tential. In addition we add a J-split in the center of the 
trap potential which may be viewed as a generic model 
for double well structure or, alternatively, as a good ap- 
proximation to the problem of a trap with an impurity 
at the center. While the the optical box trap has been 
experimentally achieved [20] , the 5-split potential may be 
realized by adding an additional laser in the center of the 
box trap. By changing the strength of the laser, one can 
tune the height of the split barrier effectively. In order 
to get some exact results in the strongly coupling limit, 
we first investigate the system by Fermi-Bose mapping 
0, ^] in the infinite repulsive limit in which the inter- 
actions between bosonic particles dominate the physics 
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of the system. Then we use the exact diagonahzation 
method to deal with the situation in the full interaction 
regime. In particular, the reduced single particle density, 
momentum distribution and occupation of the natural or- 
bits are illustrated by varying the parameters such as the 
number of particles, the strength of interaction between 
particles and the height of split barrier. 

This paper is organized as follows. In Sec. II we de- 
scribe the Hamiltonian of the few-boson system and re- 
view the single particle eigenstates. In Sec. Ill, the TG 
gas is studied and we will see how the properties of the 
ground state in this extreme limit reveal the implication 
of this few-boson system. Subsequently Sec. IV is devoted 
to the exact diaganolization method to deal with the vari- 
able interaction case. Finally we summarize our result in 
Sec.V. 



II. MODEL HAMILTONIAN AND 
SINGLE-PARTICLE EIGENSTATES 



Here, C is the normalization constant and p is the wave 
vector of particle determined by p//^ + tan {pa) = 0. The 
corresponding energy is E = p'^ /2. By contrast, the an- 
tisymmetric eigenfunctions are unaffected by the barrier 
and vanish at the origin. They are the same as the odd 
eigenstate of the hard wall trap without split {hz = 0) 
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The corresponding energies is determined by En = 
((n + l)7r/2a)^ /2. In the limit k ^ 0, energies 
for symmetric states in unit of 7r^/2a^ are E^ = 
(1/2) ,(3/2) ,(5/2) , and those for the antisym- 

metric eigenstates are = 1, 2^, 3^, • • • . In the opposite 
limit, i.e. ^ oo^ the former energy series converge to- 
wards En = 1, 2^, 3^, • • • , and each symmetric eigenstate 
becomes degenerate with the next highest-lying antisym- 
metric state. 



We consider TV bosonic atoms of mass m confined in a 
highly elongated hard wall trap which can be treated as 
quasi-one dimensional. For the case study of localization 
of particles, we split the trap by inserting a tunable zero- 
ranged barrier at the trap center [18]. The atoms interact 
via a delta potential and the Hamiltonian for this system 
is {h = m = 1) 
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where the single-particle Hamiltonian hi is given by 



hi 



1 

'2dxJ 



V (Xi) + K.S (Xi) . 



(2) 



Here V (x) describes a hard wall trap which is zero in the 
region (—a, +a) and infinite outside. The last term in the 
single-particle Hamiltonian represents a 5-type barrier lo- 
cated at the origin and we denote the strength of this 
barrier by a positive parameter k,. Another essential pa- 
rameter to characterize the system is the one-dimensional 
interaction strength g which is determined by the modi- 
fied s-wave scattering length due to the strong transverse 
confinement [Sj]. The single-particle Hamiltonian can be 
exactly solved and the solution can be found in standard 
textbook and we review the result here again for later 
computation. 

The eigen equation of single-particle Hamiltonian is 
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• hzd {x) 



^{x)=E^{x). (3) 



The eigenfunctions are either symmetric or antisymmet- 
ric. The analytic symmetric eigenfunctions are 

rh (^\- I " P '''' -a<x<0 



III. TONKS-GIRARDEAU GAS 

The infinite interaction between the atoms can be rep- 
resented as a constraint on the allowed bosonic wavefunc- 
tion ^(xi, • • • , xn). 



,Xn) 



if Xi = Xj,l <i < j < N, (6) 



which suggests that TG gas may be viewed as a group 

N 

of "free" bosons governed by the Hamiltonian ^ hi, 

i=l 

while its wavefunction is subjected to the constraint in 
Eq. (|6]). One can see immediately that this constraint 
is equivalent to the exclusive principle for a Fermi sys- 
tem with antisymmetric wavefunctions. Based on this 
observation Girardeau [21] found the Bose- Fermi map- 
ping theorem that allows one to construct wavefunction 
for a TG gas from the wavefunction of a free Fermi 
system. In particular, the wavefunction of the bosonic 
system ^b(^i, • • • , ^at) is related to the non-interacting 
fermionic wavefunction ^i?(xi, • • • , xat) by the mapping 



^b(^i. 



,xn) = A{xi,--- ,xn)'^f{xi, - ■ ■ ,xn)- (7) 



The fermionic wavefunction can be constructed 
from Slater determinant ^i?(xi, • • • , xat) = 
y^l/N\ det^Jo^^^i (^n {xj)) where (^j) are the 
single particle wavefunctions obtained above. The 
antisymmetric function A(xi, • • • , xat) is defined as 

A(xi, • • • , xat) = Yl^9^{^i~^j) with sgn the sign 

i>j 

function. The TG gas is exactly solvable in this way 
and one easily see that the single-particle density or the 
thermodynamic properties are identical for these two 
systems. 

However, quantum correlations considerably differ 
from each other as can be manifested in the reduced 
single-particle density (RSPD) or the momentum distri- 
bution. Although the exact many-body wave function 
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FIG. 1: (Color online) Reduced single-particle density matrix 
of 5 bosons in TG limt, for barrier strength — 0.2, {\:>)k — 
2, (c)/^ = 20, (d)/^ = 200. Each plot spans the range —a < 
x,y < a. 

describing TG gas can be written in compact form, it is 
a difficult task to calculate the RSPD defined as 

p{x^y) = N J dx2 ■ ■ ■ dxN 

x^^(x,X2,--- ,^iv)^B(^,^2,--- ,xn) (8) 
and the momentum distribution 

n{k) = / dx I dyp (x, y) exp {-ik {x - y)) 

(9) 

which is normalized to one. To resolve the difficulty 
of multi- integration, Pezer and Buljan Q rewrote the 
RSPD in terms of the single particle basis 

N 

p (x, y)=Y. ^^ (^) y) ^3 {y) ' (10) 

The coefficients Bij (x, y) constitute an A/" x A/" matrix 
B (x, y) = (P~^) det P and the entries of the matrix P 
are Pij (x, y) = 5ij — 2 Jj^ dx'cl)* (x^) {x')^ assuming x < 
y without loss of generality. This scheme enables efficient 
and exact numerical calculation of quantum correlations 
in few body systems. 

We illustrate the RSPD of 5 particles in Fig. 1 and 
that of 6 particles in Fig. 2 for four different barrier 
strengths n = 0.2, 2, 20, 200, respectively. The diagonal 
term p (x, x) is nothing but the single-particle density 
and is normalized to N . For a negligibly small split bar- 
rier, Fig.l (a) and Fig. 2 (a) show that the diagonal line 
have N peaks meaning N density maxima, since the in- 
finite interaction between particles repel each other to 
occupy different positions. Meanwhile the off-diagonal 
elements are non-zero, which reflect the correlation be- 
tween particles. Physically this self correlation p{x^y) 
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FIG. 2: (Color online) Reduced single-particle density matrix 
of 6 bosons in TG limt, for barrier strength {a)K = 0.2, {h)K, = 
2, (c)Hi = 20, (d)Hi = 200. Each plot spans the range —a < 
x,y < a. 



can be viewed as the probability that two successive mea- 
surements will find the particles at position x and re- 
spectively. Increasing the barrier strength, as seen in 
Fig.l (b-d) and Fig. 2 (b-d), leads to the emergence of a 
quadrant separation which coincides with the result of 
Ref. [18|]. We note also in this system there appears an 
interesting parity effect for odd and even number of par- 
ticles. When N is odd, the middle peak begins to split 
into two with increasing barrier strength, so the number 
of density peaks in each well increases from {N — 1) /2 
to {N + 1) /2. When TV is even, on the other hand, there 
are always N/2 peaks in each well for any value of barrier 
height. For a very strong split in the center, e.g. hi = 200, 
the off-diagonal quadrants diminish for N is even, while 
remain nonzero for N is odd. This implies that for odd 
number of particles, the correlation between the two sides 
of the split persists for a rather high barrier. 

We can explain our result in a way similar to the energy 
band theory. The barrier in our system splits the hard 
wall trap into a double- well structure. Just as the exten- 
sion of two-site Hubbard model to a lattice gives natu- 
rally the energy band, one can imagine that the hard wall 
trap with periodically distributed J— splits, i.e., Kronig- 
Penny potential, admits also similar band structure [23]. 
Each band in our system, however, consists of only two 
energy levels. The infinite interaction between particles 
forbids particles to occupy the same level and only two 
particles are allowed to populate in each energy band. 
For the ground state, particles occupy the energy lev- 
els from the bottom one by one. When N is even, all 
occupied energy bands are full and the system is in an 
insulator state, while for N is odd, the highest occu- 
pied band is half filled which reminds us the conduction 
band in metal conductor. Particles occupying that band 
can move more freely which is the reason for different 
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FIG. 3: (Color online) Momentum distribution, n{k), of 5 par- 
ticles (a) and 6 particles (b) in TG limit for k = 0.2, 2, 20, 200. 
When the number of particles is odd, the secondary peaks ap- 
pear for increased barrier strength due to the interference of 
particle in two almost separate wells, while for even, the mo- 
mentum distribution broaden for increased barrier strength 
owing to the localization of particles in separate halves of the 
trap. 
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FIG. 4: (Color online) Occupation of the natural orbits Xi 
of 5 particles (a) and 6 particles (b) in TG limit for k = 
0.2, 2, 20, 200. When the number of particles is odd, with the 
increased barrier strength, the two largest occupations do not 
change almost while the other nearby two occupations have 
the same value which form pairs. When N is even, with the 
increased barrier strength, all nearby two occupations have 
the same value which also form pairs. 



oflF-diagonal elements of RSPD induced by the parity of 
particle number. 

These results can also be understood in a simpler way. 
When N is even, each well can keep N/2 particles, and 
the chemical potential can keep balance in two sides of 
middle barrier. When N is odd, supposed each well is 
firstly filled with {N — 1) /2 particles, one would face a 
dilemma on how to put the last one. The probability of 
finding the last particle is the same in each well, which 
keeps the connection between the two wells and leads to 
non-zero RSPD in off-diagonal quadrant. 

To characterize some BEC-like coherence effects in few 
body system, one usually examine the eigenvalues and 
eigenfunctions of the RSPD [22[. The eigenfunctions of 
the RSPD, called the natural orbits cp {x) in TG systems 

J dyp (x, y) (fi (y) = Xi^pi (x) , z = 0, 1, 2 • • • (11) 

represent a set of effective single-particle states, while 
eigenvalues represent their occupation and = TV. 

The momentum distribution is related to these natural 
orbits via the relation n{k) = Xi\fii (k) p, where fii {k) 
denotes the Fourier transformation of cpi (x) 

jd.^.iy)e.H-^k.)/V^. 

Fig. 3 shows the momentum distribution in the TG 
limit for four different values of the barrier strength 
n = 0.2,2,20,200. When N is odd, as the barrier 
strength is increased, the main peak firstly increases, 
reaches the critical value for a barrier strength, and then 



begins to decrease thereafter. The secondary peaks ap- 
pear with increasing strength of barrier. The barrier sep- 
arates the particles into two half wells, which, however, 
are not independent but interfere with each other. It is 
this interference that leads to the appearance of the sec- 
ondary peaks. When the number of particles increase, 
the secondary peaks become more and more slight as a 
result of infinite interaction between particles. When N 
is even, the momentum distribution becomes broader and 
broader with increasing barrier strength and the height of 
peak decreases monotonically, which coincides with ear- 
lier observation that N/2 particles become localized in 
each separate half-well and keep balance in the two sides 
of barrier. 

We also present the occupation of the natural orbits 
for odd and even particles in Fig. 4. Clearly the presence 
of the split in the center will not affect the occupation 
obviously for odd number of particles. The first and sec- 
ond occupancies remain the same for increasing barrier 
height, but the third and fourth occupancies tend to form 
the staircase structure, implying the formation of pairs. 
Subsequently, the next two occupancies also form pairs, 
and so on. Contrary to this, for even number of parti- 
cles, all occupations tend to form pairs with increasing 
barrier strength due to the degeneracy of energy levels, 
which is also observed in fermionic TG gas [25]. In the 
case of odd number of particles, the presence of the last 
particle prevents the formation of pairs in the the first 
and second occupations. 
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FIG. 5: (Color online) Particle density of five bosons for 
barrier strength hi — 2 and different interaction strength 
^ = 0, 2, 5, 20, 100. The density becomes flatter and broader 
with increasing g. In the strong interaction regime new peaks 
appear gradually and finally the density tends to be the same 
as that of 5 free Fermions. 

IV. VARIABLE INTERACTION STRENGTH 

For finite particle interaction, we use exact diagonal- 
ization method [14] to deal with the many body Hamil- 
tonian ([T|), which takes the second quantized form 

^ = ^ Eia\ai + ^ Iijkia\a]akai (12) 

i ijkl 

where a] (a^) is the bosonic creation (annihilation) oper- 
ator for a particle in the single particle energy eigenstate 
(pi {x). The interaction integration parameters lijki are 
calculated through lijki = / dx(j)i {x) {x) (j)k (x) 4>i (x). 
The Hamiltonian is diagonalized in the subspace of the 
energetically lowest eigenstates of non-interaction many- 
particle system. Some exact results for a pair of bosons 
have been obtained in ref. [18] and our calcualtion for 
even number of particles implies similar behaviours in 
the relevant physical quantities. We will concentrate on 
the results for odd number of bosons = 5 considering 
26 single particle eigenstates in the following. 

Let's first discuss the particle density, p (x) = 
<^^t(^)^(^)^ = Ei,j (^)^i (^) where the 

mean value is calculated on the ground state of Hamil- 
tonian (p!2|) . Without loss of generality, we choose the 
barrier strength k = 2. Fig. 5 shows the particle density 
for variable interaction strength ^ = 0, 2, 5, 20, 100. Be- 
cause of the presence of the central barrier, the density in 
the middle position appears to be lower than its nearby 
region. For noninteracting bosons, there is one peak on 
either side of the barrier respectively. As the interac- 
tion strength is increased, the particle density begins to 




FIG. 6: (Color online) Momentum distribution of five bosons 
for different interaction strength g. The secondary peaks dis- 
appear with increasing interaction strength which hinders the 
interference of the particles in two separate half-wells. 



spread out and becomes flatter, while the height of two 
peaks is decreased. New small peaks appear for certain 
interaction and finally the density tends to be the same 
as that for free fermion system. 

The momentum distribution for finite interaction is 
now calculated as n{k) = (fl'^ (k) fl (/c)^ /N and is nor- 
malized to one. Here H (k) annihilates a boson with 
momentum k and can be expanded upon the Fourier 
transformation of single particle eigenstates {x). Fig. 6 
shows the momentum distribution for the same values of 
interaction strength as in Fig. 5. We find, as the inter- 
action is increased, the secondary peaks disappear grad- 
ually and the central peak becomes lower and broader. 
This result proves that the interaction will diminish the 
interference between two wells. The momentum distribu- 
tion is identical for ^ = 100 and TG limit g = oo which is 
calculated by the Fermi-Bose mapping in the above sec- 
tion. In the TG case, the long tail appears because the 
higher energy levels may been occupied by some particles 
owing to the infinity repel interaction, which agrees with 
general nature of a TG gas. This is in strike contrast 
with the main feature of momentum distribution for 5 
free fermions which takes the plateau shape composed of 
5 peaks and without the long tail. 

We also present the result for occupation number dis- 
tribution {rii) = (a'lai) which is different from the occu- 
pation of natural orbits mentioned above. The occupied 
states here are single particle wavefunctions in the trap, 
while natural orbits are eigenfunctions of the RSPDM. 
As can be seen in Fig. 7, with increasing g, the bosons 
occupy not only the ground state but also some low-lying 
excited states. The occupation number in the lower state 
are always larger than that in the next higher one, which 
is different with particles in the harmonic trap [l3]. In 
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FIG. 7: (Color online) Occupation number distribution of five 
bosons for different interaction strength g. On each state in 
(d), the left solid part represents g = 100 and the right hollow 
part shows the TG case. With increasing g the bosons leave 
the ground state and occupy excited states. The occupation 
number in the lower state are always larger than that in the 
next higher state. 



FIG. 9: (Color online) Momentum distribution of five bosons 
for varying inter-particles interaction strength g = (a), 
2 (b), 20 (c), oo (d). Within each plot the distribution 
is considered for 4 different values of the barrier strength: 
K = 0.2,2,20,200. Secondary peaks appear with increasing 
barrier strength, while disappear for increasing interaction 
strength. 




FIG. 8: (Color online) Particle density distribution of five 
bosons for varying inter-particles interaction strength g = 
(a), 2 (b), 20 (c), oo (d). Particles are separated into two 
parts with increasing barrier strength. 



the harmonic trap, single-particle states with even parity 
are stronger occupied compared to those with odd par- 
ity. The difference between our double well and single 
well structure comes from the fact that odd parity states 
can never feel the existence of the split in the center. 
In Fig. 7 (d) we compare these two kinds of occupations 
of single particle states and natural orbits and find that 
they almost have the same value in the corresponding 
state though the details are not completely identical. 



The last knob to control the system is the variable 
barrier strength of the split. Fig. 8 shows the parti- 
cle density distribution for four inter-particle interac- 
tion strengths ^ = 0, 2, 20, oo and four values of barrier 
strengths = 0.2, 2, 20, 200. As the barrier strength is 
increased, the density in middle area is decreasing while 
the density in other area is increasing for any interac- 
tion. The barrier therefore divides the particles into two 
areas regardless of the interaction. It is more construc- 
tive to say more words on the secondary peaks in the 
momentum distribution, which has already been shown 
in Fig. 6. Physically, the interference between particles 
split into the two separate half-wells gives rise to these 
peaks, which reminds us the famous double-slit arrange- 
ment. Fig. 9 (a) shows the momentum distribution of 5 
non-interacting particles {g = 0) for varying k,. In this 
case the momentum distribution is given by the square 
of the single-particle wave-function in momentum space. 
So we can write it out analytically 



n(p) = C{ 



sin {k + p) sin {k — p) 



k -\- p 

hzl — COS {k 

k 



k 



hi 1 
k~ 



-P) 



p 



cos {k + p) 



(13) 



with C is a normalization factor. The secondary peaks 
can be seen more and more clearly with increasing 
strength of barrier. Strength of the barrier causes the 
spreading of the particle-density in two wells, which in 
turn leads to a narrowing of momentum peaks. For finite 
interaction strength. Fig. 9 (b-d) show similar tendency 
for the central and secondary peaks. The disappear- 
ance of the secondary peaks with increasing interaction 
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strength can be attributed to the increased locahzation of 
the particles, analogous to the superfluid-insulator transi- 
tion of the Bose-Hubbard model in optical lattice studies 




V. CONCLUSION 

In conclusion, we have studied the ground state prop- 
erty of few bosons in a hard wall trap with (5— type split 
in the center. Based on the Fermi-Bose mapping, the 
reduced single-particle density and the momentum dis- 
tribution of odd {N = 5) and even {N = 6) bosons are 
illustrated and explained in TG limit. The off-diagonal 
quadrants representing the interference between two half 
wells remain nonzero for N is odd, while disappear for N 
is even. For finite interaction strength, the exact diago- 



nalization method enables us to manipulate the system 
in several ways, including variable number of particles, 
interaction strength, barrier height of the split, etc. The 
interference between particles on the two sides of the split 
manifests itself in the secondary peaks of the momentum 
distribution, emerging with increasing barrier strength hz, 
while vanishing for increasing interaction strength g. 
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